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Abstract 

We analyze the service times of customers in a stable M/M/l queue in equilibrium 
depending on their position in a busy period. We give the law of the service of a 
customer at the beginning, at the end, or in the middle of the busy period. It enables as 
a by-product to prove that the process of instants of beginning of services is not Poisson. 
We then proceed to a more precise analysis. We consider a family of polynomial 
generating series associated with Dyck paths of length 2n and we show that they 
provide the correlation function of the successive services in a busy period with n + 1 
customers. 

Resume 

On s'interesse a l'analyse des temps de service des clients d'une file M/M/l stable 
et en equilibre selon leur position dans une periode d'activite. On donne la loi d'un 
service sachant que le client se trouve au debut, a la fin ou au milieu de la periode 
d'activite. Ceci permet, au passage, de prouver que le processus des instants de debut 
de service n'est pas un processus de Poisson. On mene ensuite une etude plus fine. On 
exhibe une famille de series generatrices polynomials associees aux chemins de Dyck 
de longueur 2n et on montre qu'il s'agit de la fonction de correlation des diffcrents 
services dans une periode d'activite comportant n + 1 clients. 

Keywords: M/M/l queue, busy period, Dyck paths. 
AMS classification (2000): 60K25, 68R05. 

1 Introduction 

The M/M/l/oo/FIFO queue (or M/M/l queue) is the queue with a Poissonian arrival 
stream, exponential services, a single server, an unlimited buffer capacity, and a First-In- 
First-Out service discipline. It can be argued that the M/M/l queue is the most elementary 
and the most studied system in queueing theory, see for instance O [HI El E] . Quoting [9] , 
"most likely, any book with queueing in the title has something to say on the subject". 

Let A be the intensity of the Poisson arrival process and let \x be the parameter of the 
exponential service times. Assume that the stability condition A < /i holds and consider 
the queue in equilibrium. Our objective is to get precise information on the distribution 
of the service of a customer based on its position in the busy period. 
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First of all, recall that the distribution of the first, respectively last, service is an expo- 
nential of parameter fj,, respectively (A + fj,). We are then able to compute the distribution 
of a service in the "middle" of a busy period (i.e. neither at the beginning nor at the 
end). As a by-product, we also get the distribution of the duration between two successive 
beginning of services. Since it is not an exponential, we conclude that the point process 
of the instants of beginning of services is not Poisson (as opposed to the point process of 
completion of services) . 

Then we study the service time of the k-th customer in a busy period of length n + 1 
(i.e. containing n + 1 customers). Consider a busy period conditionned to be of length 
n + 1, and let ((/i)o,...,2n+2 be the corresponding embedded queue-length excursion. Its 
trajectories are equiprobable and it is easy to see that they are in bijection with Dyck 
paths of length 2n. If we condition ((?j)o,...,2n+2 to be associated with a given Dyck path 
7r of length 2n then we observe that the law of the service time of the k-th customer is 
equal to the convolution product of + 1 exponentials of parameter A + /i where Ik is the 
length of the intersection of it with the line y = x — 2k. By summing over Dyck paths of 
length 2n, we get an expression for the joint law of the services in a busy period of length 
n+ 1. Then using elementary properties of Dyck paths, we obtain results on services within 
a busy period somewhat difficult to obtain by direct probabilistic arguments (Section 4). 
The correlation function of the services is a natural generating polynomial of Dyck paths 
following a simple integral recursion (Section 5). 

Using the combinatorial properties of lattice paths to study the busy period of simple 
queues is classical, see [H El [12] and references therein. In these articles, quantities such 
that the area swept by the queue-length process during a busy period are studied, with a 
much more involved combinatorial analysis than what is presented below for the sequence 
of services within a busy period. This should come as no surprise. The area and related 
quantities, are derived by counting in a Dyck path the number of ascents and descents of 
a given vertical coordinate (Dyck paths are lattice paths in N x N, see Section 3). On the 
contrary, the sequence of services is derived by counting in a Dyck path the number of 
ascents of a given horizontal coordinate (roughly speaking). This is in essence like working 
with generating polynomials of Dyck paths in non-commuting variables. It is therefore 
hopeless to get as precise information. 

2 In the Middle of the Busy Period 

Given a positive real random variable X with law fi, denote its Laplace transform by 
^ x (s) = # M (s) = / exp(-sx)dfi(x),s G C,Re(s) > 0. We write L[X | A] = F{X G • | A} 
for the conditional law of X given an event A. The corresponding Laplace transform is 
denoted ^x\a( s )- The convolution product of two probability distributions [i and v is 
denoted by fi * v. The indicator function of a subset A of a set is denoted by Ja- It is 
convenient to denote by Exp(a) the exponential distribution of parameter a defined by 
Exp(a)[x, +oo) = exp(— ax),x G R+. Recall that ^Exp(a)( s ) = a/(a + s). 

We consider an M/M/l queue with the following notations. Let (A n ) n ^z be the arrival 
Poisson process of intensity A. Let (T n ) n g^ be the inter-arrival times, with T n — A n +i — A n . 
Denote by (<r n ) ne z the service times of the customers. The sequence (cr n ) n is i.i.d. and 
do ~ Exp(^). We assume that the stability condition A < \i is satisfied, and we consider 
the queue in equilibrium. Let (QtjteR be the queue-length process, where Qt is the number 
of customers either in service or in the buffer at time t. 
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The state of the server can be described as an alternating sequence of idle and busy 
periods. A busy period is a maximal period during which Qt > 0. An idle period is a 
maximal period during which Qt = 0. The length \B\ of a busy period B (not to be confused 
with its duration) is the number of customers served during the busy period. Throughout, 
when we consider a generic busy period B, we denote for simplicity by ctq, . . . , o\b\-i 
and To, ... , Ti_B|_i respectively the service times and the inter-arrival times of the different 
customers in the busy period. 

Lemma 2.1. Let A n be the event that a generic busy period consists of n + 1 customers, 
then 

i i n n 

An = £r 3 <^2a j ,i = {0,...,n-l}; > X> > " W 

j=0 j=0 j=0 j=0 

The justification is easy. 

The durations of successive busy periods and idle periods are independent random 
variables. The duration of an idle period is clearly distributed as Exp (A). The distribution 
of a busy period is more complex. The next results can be found for instance in [3l 
Chapter II. 2. 2] or [111 Chapter 1.2]. The probability that a busy period B consists of 
(n + 1) customers is given by 

\n..n+l 

F{\B\=n + l} = C n {x+ ^ )2n+1 , (2.2) 

where C n is the n-th Catalan number, see $3l Let 5 n be the conditional law of the duration 
of a busy period, given that the length of the busy period is (n + 1). The Laplace transform 
of 5 n is given by 

(A + ^+i 

*** ( ' )= (A + ^ + ^+i - (2 ' 3) 

Hence, 5 n is the distribution of the sum of (2n + 1) i.i.d. r.v.'s of law Exp(A + //). 

Given two independent random variables X ~ Exp(a) and Y ~ Exp(/3), where a, (3 G 
R?j_, recall that 

L[X \X>Y] = Exp(a + /3)*Exp(a), L[X \ X < Y] = Exp(a + 0) . (2.4) 

Using elementary arguments based on the memoryless property of the exponential distri- 
bution, we get: 

do ~ Exp(/i), <7|b|_i ~ Exp(A + fj,) . (2.5) 
Furthermore, remarking that {\B\ = 1} = {o"o < ro} and using (|2.4f) . it follows that: 

£[<7 | \B\ = 1] =Exp(A + ^), £[<7 I \B\ > 1] = Exp(A + * Exp(^) . (2.6) 

Our goal is now to derive the law of a service in the middle of B, i.e. of a service which is 
neither the first nor the last one (assuming that \B\ > 2). 

Let <7* be the service of a generic customer numbered * and let B be the busy period 
it belongs to. Define the events 

(£ = {* is the only customer of B} = {\B\ = 1} 

l£f = {* is the first customer of B and \B\ > 1} 

(£i = {* is the last customer of B and \B\ > 1} 

= {* is in the middle of B and \B\ > 2} . 
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Clearly the four events are disjoint and P{6 U 2f U €i U 6 m } = 1. Since the lengths of 
successive busy periods are i.i.d., we obtain immediately that 

P{£ f } = P{£i} 
P{£ f } = P{|B| > 1}/E[|B|] 
P{6 D U6 f } = l/E[\B\}. 

Now using §F1§, we get P{|£| > 1} = A/(A + fi) and E[\B\] = [i/(fi - A). It follows that 

Clearly, oq ~ £[<7* | G; Q U 6f] and cibi-i ~ <C[c* | 6J. We deduce that 

= P{C U(£ f }* <T ,|^ U(!f ( a )+P{(£ 1 }* <T ,|^( S )+P{6 m }* <T ,| (Sni ( S ) 
= P{<2 o U<£ f }tf CT0 ( S ) + P{6 1 }tf (T|B| _ 1 ( s ) +P{£ m }*^| £m ( S ) . 



That is, 



M p-\ p \{p-\) p + \ 2A 2 



/i + s /i fi + s fj,(fi + A) /i + A + s /i(/x + A) 

After simplification of the above expression, we obtain the Laplace transform of the con- 
ditional law of o"* on the event 6 m : 

xr, t \ (2/i + a)(/z + A) 

^' £ - (s) = 2(^ + S ) (/ , + A + 5 ) - (2 ' 8) 

As a by-product, we can prove that the process of instants of beginning of services is 
not a Poisson process, in contrast with the process of completion of services (departure 
instants) which is Poisson of intensity A according to Burke Theorem j2j [7] . Let us detail 
the argument. Let be the difference between the instants of beginning of services of two 
generic successive customers numbered * and (* + 1). Let B be the busy period of * and let 
v be the first idle period following B {y ~ Exp(A)). Using (|2.5p . (|2.6p . we get immediately 
that 

£*le = K + ^lCo => I £o] = Exp(A + /x)*Exp(A) 

£*1<Sj = folgf => I <£f] = Exp(A + /i) *Exp(/x) 

£*l<s, = (<T|B|-i + f)l^ =► M£*|6i] = Exp(A + /i)*Exp(A) , 

and 

| 6m] = I 6 m ] . 

Since we have just computed ^ at \£ m (s), we deduce the Laplace transform of £*: 

(s) = A(M 2 (A + M) + M2M + A) a + A^) 
5 /j(A + a)(/i + s)(A + p + s) v ; 

We check on this expression that E[£*] = 1/A and we have E[£ 2 ] = a3 -^+^+^ a . In 
particular, we have E[£ 2 ] < E[d 2 ] = 2/A 2 , where is a generic inter-departure time. 
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3 Dyck Paths 



The Catalan numbers (C n ) nS N are defined by 



* = ^M=r^f*|. (3-D 



2n + 1 \ n J n + 1 \ n 
The generating function of these numbers is given by 



Eg 



n _ 1 - VI - 4a; 



2x 

n=0 

The first Catalan numbers are C = 1, Ci = 1, C 2 = 2, C 3 = 5, C 4 = 14, C 5 = 42, C 6 = 
132, C7 = 429, • • • . They appear in many combinatorial contexts see for instance I10j . 
In particular, C n is the number of Dyck paths of length 2n. A Dyck path of length 2n is a 
path in the lattice N x N which begins at the origin (0, 0) ends at (0, 2n) and with steps of 
type (1,1) or (1,-1). Denote by D n the set of Dyck paths of length 2n, observe that Do 
is a singleton whose element is the unique Dyck path of length 0. 

We now define a family of polynomials related to Dyck paths. Let ir £ D n and let jj be 
the line y = x — 2j, for j £ {0, ... ,n— 1} and denote by ay the length of the intersection 
of 7j with 7r (equivalently ctj + 1 is the number of lattice points common to it and 77). We 
introduce two polynomials P w and R n defined by 

re— 1 ai re-1 

Pw(yo,yi, ■ ■ ■ ,2/n-l) = \\ , RniVQiVl, ■ ■ ■ ,Vn-l) = vT ■ ( 3 - 2 ) 
i=0 i=0 

Let (-Pn)ngN and (i? n ) ne N be the two families of polynomials defined by P n = X^eD ^1 - 
^ 7rgI)n R-K. Clearly P n and Rn are homogeneous polynomials of degree n over the n vari- 
ables y ,yi, • • • ,2/n-i- 



4 The Law of the Services in a Busy Period 

Recall that the queue-length process (Qt)teM i s a continuous time Birth-and-Death process 
on N with generator P such that P n ,n+i = A, n > 0; P n ,n-i = w > 1; -P n ,m = 0, \n — m\ > 
2. Let {q n ) n £i denote the Markov chain embedded at its jump instants. More precisely, let 
T be the point process obtained as the superposition of the arrival and departure processes 
and let (T n ) n <=z be its points with the convention To = Aq. Then we set q n = Q T -. The 
transition matrix of (q n )n& is given by 

Po,i = i; Pi,i-i = , 7 ; Pi,i+i = -r~r — ' * - ^ 

A + /i A + /i 

and pj j = otherwise. 

A busy period corresponds to an excursion of (q n )nez from to its first return to 0. 
With the same numbering convention as in Section 2, the generic busy period B consists 
of n + 1 customers if and only if 

q = 0, ft > 0, i € {1, • ■ • , 2n + 1}, q 2n +2 = . (4.2) 

On this event, the (random) path with successive edges (i — 1,^ — 1), i £ {1, • • • , 2n + 1} 
is a (random) Dyck path of length 2n. We call it the Dyck path associated with B (see 
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Figure [T]). On the event {\B\ = n + 1}, all Dyck paths appear with the same probability 
(the probability of a given trajectory (q n ) n depends only on the number of increasing and 
decreasing jumps, see (14. ip ). On the event that \B\ = n + 1 and that the associated Dyck 
path is 7T G T> n , the power of y%-\ in P n is the number of customers which join the system 
between the the i-th and the (i + l)-th departures. Combining these observations with the 
fact that the time between successive transitions of (Qt)t^R axe independent r.v.'s of law 
Exp(A + n) as long as the queue is non empty, we get : 

Theorem 4.1. Given that the length of the busy period is n + 1, the conditional density of 
the random vector (o"o, . . . , o~ n ) representing the service times of the successive customers 
is 



D(y , 



,Vn) 



(A + M) 



2n+l 



r (^)(*/o+-+ 2 /„- 1 )p n(y0) . . . ,y n ^)e-( X+ ^ 



where P n is the Dyck polynomial of degree n defined in Section 3. 



(4.3) 



Q, 





n-l 



Figure 1: Dyck path associated with a busy period. 
A direct computation of the Laplace transform leads to the following : 

Corollary 4.2. Consider a random vector (do, . . . , <f n ) ~ £[(<to, . . . , o~ n ) \ \B\ 
Laplace transform is given by 



^(a ,...,a n )(so,---,s n ) 



E[f[ e~ s ^} 
i=0 
j n— 1 

— (JJ Zi)R n (z , . . .,Zn-l)Zn , 



1=0 



where z,- t = x+^+ s i ^ {0> • • • > n }> an d Rn ^ s defined in Section 3. 



n + 1). Its 

(4.4) 
(4.5) 
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Let us paraphrase the above results in a somewhat more intuitive way. In a busy period 
of length (n+ 1), the conditional law of (o~o> . . . o~ n ) is the same as the law of (5o, . . . <r n ) 
that we now describe. The law of a n is an Exp(A + (jl) independent of (do, . . . , cr n ~i). Let 
II be a r.v. uniformly distributed over T> n = {Vi, . . . ,^c n }- Conditionally on {II = 7Tj}, 
the r.v.'s aj are independent and distributed as the sum of fc*- random variables of law 
Exp(A + fj,), where (k\ — 1) is the exponent of yj in P ni . This is illustrated in Table 1. 

We now exploit the correspondance with Dyck paths. 

Let T)< 1> be the set of Dyck paths of length 2n where the first return to the axis 
{(n, 0), n G N}, after the origin (0,0), occurs at the point (2i, 0), i G {1, . . . , n}. Clearly, 
the sets D^ l> are disjoint and D n = Uf =1 T)< l> . Furthermore 

V< 1> ~ IV i x D n -i . (4.6) 

A consequence of the above is the very classical identity on Catalan numbers : 

C n = ^ ] CiCj . 

i+j=n—l 

Let R n be defined by D?„(2 0) . . . , z„-i) = R n (z , . . . , z n -\) x zq . . . z n _i, going back to 
Corollary 14.21 we have 

^(ao,...,a„){so, - ■ ■ ,S n ) = —!R n (Z0,...,Z n -l).Z n . 

where z-i = x +t+ s - ; Vi G {0, . . . , n}. We also define = 3? n (l, . . . , 1, Zj, 1, . . . , 1), then 

Proposition 4.3. On the event {\B\ = n + 1}, we have 

4 0) (z ) = z ,Vn>0, ^ 0) (z ) = ^ ^i 0) («b)Ci 5 ( 4 - 7 ) 

i+j=n— 1 

/or, < k < n — 1, 

fc-i 

a?°(*k) = C*** , %^(z k ) = Y, ^(z^Cj + E^Vfc(^) ■ (4.8) 

i+j=n— k— 1 Z=0 

On Table 1, one notices a simple relation between the laws of ctq and o~i, which is actually 
always true : 

Proposition 4.4. Let B be a generic busy period, for n > 1 we have 

£[cr | \B\ =n+l] = £[<ti | |S| = n+l]*Exp(A + ^) . (4.9) 
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Figure 2: The mapping ^ : T> n — ► D 
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28 
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Table 1. Services. 



The table should be read as follows. For instance, on {\B\ = 5}, the law of 02 is 
fi a3 = ^Exp(A + /u) + ^Exp(A + ^)*Exp(A + /i) + ^Exp(A + ^)*Exp(A + /i)*Exp(A + /x). 

Proof. The mapping \& : D n — > D„ is defined in Figure [2] It is clearly an involu- 
tion, hence a bijection. More formally, given a Dyck path tt £ T> n such that R w = 
yi-yoQ(yo,y2,---,y n -i) then ^(vr) G D n is defined by i?^ (7r) = y^ +1 Q(yi,y 2 , . . . ,y n -i)- 
In view of Corollary 14.21 it completes the proof. □ 

5 Dyck Paths Polynomials 

We go back to the family of polynomials (P n ) ng N defined in Section [3l We are going to 
use Theorem 14.11 to give nice expressions for the P n 's. Let A n be the event that a generic 
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busy period consists of n + 1 customers. Let Aq, . . . , A n be borelians of R+, 

»{(Ti G A, i = 0,...,n;A n } 



P{o-i GAi, « = 0,...,n| A„} 



P{A„,} 



Let L n = P{crj G Ai, i G {0,ra};A n } and for k = 0, . . . , n, let Yfe = X^=o2/« anc ^ -^fc 
X)i=o x « • Using Lemma 12. 1| we have 

/• /-Yo rYi— Xq 

L n = (Xfi) n+1 e-^dy ...dy n e- Xx °dx e^dxf- 

J Aox—xA„ Jo JO 



Y n -i—X n -2 poo 



' Aox—xA 

... / e'^^dxn-i I e~ AXn dx n 

JO JYn-Xn-! 
f [Yo pYl—Xo l-Yn-l—Xn-2 

= A>" +1 / e -(M+A)y„ dyo _ _ dyn / dxo / dsi . . . / dXn l 

Ja x—xA„ Jo Jo Jo 

Then, using theorem 14. 1\ we get 

ryo rvo+vi-xa ryo-\ — bi/n-i— tecH — 
P n (yo, • • • ,2/n-i) = / dx dxi... dx n -i. (5.1) 

io Jo Jo 

Simple manipulations of formula (|5.1|) then yield : 

Lemma 5.1. The polynomials (P n ) n m satisfy the following equations 

/■ 2/0+2/1 

P n (yo,---,yn-i)= Pn-i(y,y2,---,y n -i)dy (5.2) 

J yi 

and 

<"yo ryi+%0 ry n -i+x n -2 



P n (yo, ■ ■ ■ ,2/n-l) = / dx dxx ... dx n ^i . (5.3) 

Jo Jo Jo 

For completeness, here is a direct proof of (|5.2|) without using Theorem 14.11 

Let Dn be the set of all Dyck paths of length In starting with i steps of type (1, 1) 

followed by one step of type (1, —1) and define the polynomial Pn^ such that 



Clearly, we have 



Pn(y , yn-l) = E fPnHvi, ■ ■ ■ , 0n-l) • 

Hence, we get 
ryo+yi ryo+yi y« ^ 

/ Pn-x(y,y2,---,y n -i)dy = / 2j-T p n-i(2/2,---,yn-i) ^2/ 

Jyi Jyi i=i l ' 

n-l 



n-l 



<=1 ^ + J - 

n 7 n— j t> 

1^7\l^~k\ p n-i m,---,yn-i) 

j=l J ' k=0 
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V 



i i \M VN SA > 



2(n-l) 



II 



2n 



Figure 3: Proof of the equality (|5.4p : Paths contributing to P^ n (Left) and to P^ 
(Right). 

With the help of Figure [3l we notice that 



n-l fe-i+1 n-3 k 



fe=j- 



■ x (k-j + 1)1 



k=0 



(5.4) 





Figure 4: The volumes of the gray areas are P2(yo,Ui) (left) and Ps(yo, yi, y<z) (right). 



It leads to 
rvo+m 



i=l 



Pn-i(y,V2, ■ ■ ■ ,yn-i) dy = ^2^-P$(yi, ■ ■ ■ ,y n -i) = P n (yo, ■ ■ -,y n -i) ■ 



This result can also be proved using the theory of species presented in [T]. Finally, using 
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(|5.3p . the polynomials (-P n ) n eN can be interpreted as volumes. We give a representation of 
this in Figure U] for n = 2 and n = 3. 

Conclusion. Here are several other simple models of queues for which the queue-length 
process is a Birth-and-Death process: the M/M/K/oo queue, the M/M/oo queue, or the 
M/M/K/L queue {K < L < oo). In each case, if the generic busy period is of length n + 1, 
we can associate with it a Dyck path of length In. However, the different Dyck paths of 
length 2n are not equiprobable anymore. Hence, we do not get a simple formula for the 
joint law of the services as in Theorem 14.11 
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